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In  this  paper  we  present  a  theoretical  analysis  for  the  interaction  of  a 
high-energy  laser  beam  with  a  hygroscopic  aerosol.  Droplet  absorption 
simultaneously  heats  the  atmosphere  and  vaporizes  the  droplet.  When  a  droplet 
is  vaporized  a  thin,  non-equilibrium  region  exists  at  the  droplet  surface.  This 
interfacial  discontinuity  determines  the  amount  of  droplet  superheat  as  a  function 
of  the  rate  of  vaporization.  For  submicron  droplets  the  degree  of  superheat 
achieved  is  small  because  heat  conduction  to  the  atmosphere  keeps  the  droplets 
cool.  Less  than  ten  percent  of  the  absorbed  energy  is  available  for  vaporization. 
Droplets  larger  than  l|am  are  heated  to  higher  temperatures  with  as  much  as 
ninety  percent  of  the  absorbed  energy  available  for  vaporization.  After  the 
large  droplets  are  vaporized  to  submicron  size  the  residual  dropltets  scatter, 
absorb  and  thermally  bloom  the  laser  beam  as  if  they  are  non -vaporizing. 
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I. 


INTRODUCTION 


There  have  been  several  theoretical  and  experimental  investigations  on 

the  interaction  of  a  high-energy  laser  beam  with  hygroscopic  aerosols.  This 
problem  is  of  interest  because  the  absorbed  energy  can  vaporize  the  droplets  and 
increase  the  meteorlogical  range.  However  energy  transport  between  the  heated 
droplet  and  the  atmosphere  simultaneously  occurs  with  vaporization  and  reduces 
the  energy  available  for  vaporization.  In  this  paper  we  present  the  analysis 
for  determining  the  fraction  of  absorbed  energy  that  results  in  vaporization.  The 
remainder  of  the  absorbed  energy  is  deposited  into  the  atmosphere  as  heat  and 
causes  thermal  blooming.  To  accomplish  this  we  require  a  rather  detailed  treat¬ 
ment  of  both  the  thermodynamics  and  hydrodynamics  of  a  vaporizing  droplet. 

When  a  droplet  is  vaporized  a  very  thin  but  finite  non  -  equilib  rium  region  exists 
at  the  droplet  surface  to  allow  for  the  phase  transaction.  This  interfacial  dis¬ 
continuity  requires  a  saturated  vapor  pressure  at  the  droplet  temperature  that 
is  higher  than  the  actual  vapor  pressure  at  the  droplet  surface.  The  droplet  is 
said  to  be  superheated  with  the  degree  of  superheat  dependent  upon  the  rate  of 
vaporization.  The  non- equilibrium  theory  requires  high  droplet  temperatures 
and  increases  the  energy  deposited  into  the  atmosphere  as  heat.  Energy  trans¬ 
port  from  the  droplet  to  the  atmosphere  occurs  either  gas  rlynamically  or  by 
f^iff.i cion.  Convective  or  gas  dynamic  vaporization  occurs  wiienever  the  vapor 
pressure  at  the  droplet  surface  is  grecil«.r  tl'"' n  the  anrbient  press\ire.  As  e%'ap- 
oration  proceeds  the  droplet  becomes  smaller  and  the  absorbed  ir radiance, 

drcjplct  Icmipe  rature,  and  vapor  pressure  all  decrease.  When  the  v-apor 
pressure  at  the  droplet  surface  bec(>mes  less  than  the  ambient  pressure  cliffusion 

determines  the  rate  of  energy  transport  into  the  a  tuoos  j)he  re .  For  subinicron 
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droplets  the  dej^rec  of  superheat  required  for  rapid  vaporization  cannot  be 
achieved  because  heat  conduction  into  the  atmosphere  keeps  the  droplets  cool. 


Vaporization  of  submicron  droplets  is  said  to  be  kinetically  limited  with  less  than 
ten  percent  of  the  absorbed  ene rgy  available  for  vaporization.  Vaporization  of 
large  droplets  is  not  kinetically  limited  and  up  to  ninety  percent  of  the  absorbed 
energy  is  available  for  vaporization.  Thus,  an  order  of  magnitude  larger 
irradiance  is  required  to  evaporate  a  submicron  droplet  at  the  same  relative 
rate  as  a  lOlim  droplet.  Even  if  the  energy  for  complete  vaporization  is  available 
the  non-equilibrium  theory  of  evaporation  imposes  a  limit  on  the  final  droplet 
size  that  can  be  attained.  A  high-energy  laser  beam  propagates  through  the 
resulting  haze  as  if  the  droplets  are  non- vaporizing.  The  residual  droplets 
scatter,  absorb  and  thermally  bloom  the  laser  beam. 

Because  the  droplets  are  superheated  they  can  be  shattered  rather  than 

vaporized.  The  same  nucleation  phenomena  that  occurs  in  boiling  and  cavitating 

13 

flows  can  also  occur  in  heated  droplets.  Nucleation  is  the  forniation  of  vapor 
cavities  in  the  interior  of  the  droplet.  The  droplet  will  be  shattered  when  the 
saturated  vapor  pressure  at  the  droplet  temperature  exceeds  the  limiting  tensile 
stress  of  the  droplet.  The  temperature  at  which  vapor  cavities  are  formed  is 
called  the  incipient  boiling  temperature.  At  one  atmosphere  this  temperature 
varies  from  a  few  degrees  in  excess  of  100°C  when  entrapped  air  is  present  to 
267°C  in  pure  v/ater^^. 
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II.  VAPORIZATION  MODEL 


The  interaction  of  a  laser  beam  with  aerosol  droplets  that  are  vaporized  is 

!  4  .  . 

complicated  by  non-equilibrium  thermodynamics  .  When  two  phases  exist  in 

equilibrium  the  second  law  requires  continuity  of  temperature  and  the  free  energy. 

However  with  a  finite  rate  of  mass  transport  the  two  phases  are  no  longer  in 

equilibrium  and  a  thin  transition  layer  is  required.  Experimental  evidence  indi- 

-6 

cates  that  the  thickness  of  the  layer  is  a  few  mean  free  paths(10  cm)  and  is 

accompanied  by  steep  gradients  in  the  vapor  pressure.  The  net  mass  flux  across 

the  interface  is  the  difference  between  the  Langmuir  rate  of  evaporation  of  a 

pure  liquid  into  a  vacuum  and  the  rate  of  condensation  of  vapor  molecules  on  a 

1  5 

liquid  surface  and  has  been  given  by  Schrage  as 


=  e(p^(Tp)  -  F(v)p)/(2nR^Tp)‘ 


where  e  is  the  evaporation(=  condensation)  coefficient,  p®  is  the  saturated 
vapor  pressure  of  the  droplet,  p  is  the  actual  vapor  pressure  of  the  droplet, 

Tp  is  the  temperature  of  the  d  roplet  and  =  0.4621/ gm-°K  is  the  gas  constant 
for  water  vapor.  The  function  F(v)  represents  the  influence  of  the  bulk  velocity 
moving  away  frcun  the  droplet  and  is  given  by 


F(v)  =  cxp{-v‘'/2R  T„)  -  (n/2R^T^)  v(  1  -  e  rf  (v  /  (2R  T  )^));  v  =  m/p. 

v  F  F  P 


A  small  tempo  ra  lu  1 1-  jump  also  occurs  across  the  licjuid -vapor  interface  but 


it  iias  been  neglected.  The  physical  exp!anati<m  of  Eq.  (1)  follows  from  tlic  kinetic 
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theory  of  gases.  As  the  droplet  evaporates  the  vapor  molecules  collide  with 
other  molecules  and  attain  a  rando.a  velocity  distribution.  Some  of  the  vapor 
molecules  also  collide  with  the  droplet  and  reattach.  The  fraction  of  vapor 
molecules  that  stick  to  the  droplet  is  the  condensation{  =  evaporation)  coefficient 
with  the  remaining  vapor  molecules  reflected  from  the  droplet.  The  published 
values  of  the  evaporation  coefficient  vary  from  1.0  in  a  vacuum  to  0.  04  at  one 
atmosphere  ^  When  the  bulk  velocity  equals  or  exceeds  the  mean  thermal 
speed  (ZR^Tp)^  F  is  small  because  few  vapor  molecules  have  a  velocity  in  the 
direction  opposite  to  the  bulk  velocity  that  is  large  enough  to  cause  a  collision 
and  subsequent  reattachment  to  the  droplet. 

Net  evaporation  occurs  when  the  saturated  vapor  pressure  p^(T^)  of  the 
droplet  exceeds  the  actual  vapor  pressure  in  the  gas.  Net.  condensation  occurs 
when  the  saturated  vapor  pressure  of  the  droplet  is  less  than  actual  vapor  pressure 
The  non  - equilib rium  theory  requires  superheating  of  the  vapor  for  evaporation 
and  subcooling  of  the  vapor  for  condensation.  After  replacing  p  by  p/R^Tp 
in  m  we  can  dete rmine,from  Eq.(l),  the  ratio  of  the  saturated  vapor  pressure 
p^{Tp)  to  the  actual  vapor  pressure  p  as  a  function  of  the  Mach  number 
v/(YR^Tp)^  of  the  vapor.  This  is  shown  in  Fig.  1  .  For  an  evaporation  coeffi¬ 
cient  of  1  significant  non ■  equilibrium  effects  occur  for  Mach  numbers  greater 
than  0.  5.  For  an  evaporation  coefficient  of  0.  04  non  -  equilib  ri  vim  effects  are 
dominant  for  virtually  all  Mach  numbers.  Also  shown  in  Fig.  1  are  tlie  asymp¬ 
totes  where  r ccondensation  is  negligible  and  the  vacuum  rate  of  evaporation 
applies. 

Under  the  assumption  that  the  temperature  is  uniform  across  the  droplet 
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mass  conservation  and  an  energy  balance  at  the  droplet  surface  gives  the 
following  relationships  between  the  droplet  radius,  temperature  and  rate  of 
evaporation; 


p  a  =  -m  ;  m  (L  +  h  (T  )-h  (T  )  +v  /2)+  KT^  +  op^CTjZ  =  IQ  /4 
p  vpv  rpi-  a 


(3) 


^  • 

where  a  is  the  droplet  radius,  p^=  1  gm/cm  is  the  density  of  the  droplet,  m  is 
the  rate  of  vaporization,  C  =4.  2J/gm-°K  is  the  specific  heat  of  water,  the 

droplet  temperature,  L=2500J/gm  is  the  latent  heat,  h^  is  the  enthalpy  of 
the  vapor,  v^  is  the  convective  velocity  of  the  vapor,  K  =2.  55x10  ^J/cm-sec-  K 
is  the  thermal  conductivity  of  air,  I  is  the  irradiance  and  Q  is  the  Mie  absorp- 
tion  coefficient  of  the  droplet.  When  the  latent  heat  dominates  the  other  terms 
in  Eq.  (3)  most  of  the  absorbed  irradiance  will  result  in  droplet  vaporization. 

To  determine  when  this  occurs  we  must  first  solve  the  conservation  equations 
in  the  atmosphere  for  the  vapor  pressure  p  and  temperature  gradient  T^.  This 
is  discussed  in  the  next  section. 
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III. 


MASS  AND  ENERGY  TRANSPORT 


We  require  two  separate  analyses  for  determining  the  rate  of  vaporization 
as  a  function  of  the  vapor  pressure  and  the  droplet  temperature.  A  diffusion 
analysis  that  is  valid  when  the  pressure  at  the  interface  is  less  than  one  atmos¬ 
phere  and  a  gas  dynamic  model  that  is  valid  when  the  vapor  pressure  is  equal 
to  or  greater  than  one  atmosphere.  In  the  diffusion  analysis  the  kinetic  energy 
of  the  vapor  may  be  neglected  and  in  the  gas  dynamic  analysis  heat  conduction 
may  be  neglected.  These  two  situations  are  shown  in  Fig, 2. 

Our  primary  interest  is  in  propagation  where  the  droplets  are  heated 
for  times  much  longer  than  the  diffusion  time.  Droplet  evaporation  can  then  be 
described  by  a  quasi-steady  analysis.  In  the  notation  of  Ref.  17  the  equations  for 
mass  and  energy  conservation  are 


"  ( 

/r2cD\^\  0  , 

mo^c''dl  --  i| 

f4) 

dr\ 

>  1-X  /dr  / 

P  d r  dr' 

^  dr  / 

subject  to 


X(rua)  ^  X{T  )  ,  T(r=0)  =T  .  X(r=oo)  =  X  ,  T(r^a=)  .  T  , 

P  P  CB 


(5) 


where  X,  the  mole  fraction,  is  the  ratio  of  the  partial  pressure  of  water  vapor 

3 

to  the  ambient  pressure,  c  (gm-moles/cm  )  is  the  molar  density  of  the  mixture, 

o  Vo 

D--0.  23cm‘^/s  is  the  diffusicjn  caefficient  and  =  2J/gm-  K  is  the  specific  heat  of 

* 

1  7 

water  \'apor.Tlie  solutions  to  the  governing  equa  ti  ons ,  urul  e  r  the  assuniption 
that  '<  and  cD  are  constant,  give  the  heat  and  nia  s  s  flux  as 
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KT^  =  ^CJT„-T  )/{exp(r4iC%/'<  )-l). 

1  ^  p  CO  p 


m  =  (pD/R  T  o)  ln[(l-XJ/(l-X)] 

*»  V  oo 


(^>) 

(7) 


Equating  Eq.  (7)  to  Eq.(l)  determines  X  as  an  implicit  function  of  the  ratio  of 

the  saturated  vapor  pressure  to  the  ambient  pressure.  After  approximating 

i 

F(v)  in  Eq.  (1)  by  (1  - v(tt/ ZR^T^) ^),  we  obtain 


2  i 

X^(T^)  =  X  +  [D(2nT„/R,  T  )  (1 -0.  5e)/Ga  ]  ln[(l -X  )/(l -X)  ]  (8) 

P  p  V  CO  CO 


where  the  temperature  dependence  of  X  (T  )  is  given  by  the  Clasius-Clapeyron 

P 

relation.  Numerical  results  are  generated  by  regarding  the  m^le  fraction  X  as 
the  independent  variable.  For  a  given  value  of  X,  less  than  unity,  we  determine 
the  heat  flux,  the  mass  flux,  and  the  temperature  from  Eqs.  (6),  (7),  and  (8). 
These  results,  when  substituted  into  Eq.  (3),  determine  the  absorbed  irradiance. 
Because  of  the  quasi-steady  assumption  the  droplet  heat  capacity  term,  ,  is 
set  to  zero  in  Eq.  (3).  The  amount  of  droplet  superheat  is  determined  by  the 
difference  X  (Tp)-X  in  Eq.  (8)  and  increases  as  the  mass  flux  increases  or  the 
droplet  radius  decreases.  For  droplets  greater  than  10pm  the  difference  is 
negligible  and  thermodynamic  equilibrium  applies.  For  droplets  less  than  lum 
the  difference  is  large  and  the  droplet  will  not  reach  surface  equilibrium.  As 
O— *0,  for  a  fixed  droplet  temperature,  the  mass  flux  approaches  a  value  indepen¬ 
dent  of  the  radius  whereas  the  heat  conducted  into  the  atmosphere  increases  “ 

like  1  /  o.  Thus  we  will  aUvays  reach  a  final  dropUft  size  where  virtually  all  of  the 
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absorbed  energy  heats  the  atmosphere.  This  is  in  contrast  to  the  conclusion 
reached  when  surface  equilibrium  is  assurred;  namely,  that  most  oi  the  absorbed 
energy  vaporizes  the  droplet.  The  importance  of  non- equilibrium  effects  has 
been  emphasized  in  Refs.  1  and  5  and  the  above  steady  state  diffusion  analysis 
is  due  to  Williams^. 

As  the  partial  vapor  pressure  XCT^)  approaches  unity  the  diffusion  model 
is  no  longer  valid  and  must  be  replaced  by  the  convective  model  shown  in  Fig.  2. 
Initially  the  vapor  pressure  exceeds  the  ambient  pressure  and  drives  a  shock 
wave  in  to  the  atmosphere.  The  discontinuity  occuring  across  the  shock  w'ave 
relates  the  pressure  at  the  droplet  surface  to  the  vajjor  Mach  number 

1  1 

M^.  =  0.  517(373/Tp)%/p^-l)/[l  +  0.  857(p/p„-l)]^.  (9) 

As  time  increases  the  vapor  bubble  expands  more  slowly  than.  *^hc  shock  and 
remains  seperated  from  the  shocked  air  by  a  contact  discontinuity.  When  R^>a 
the  infhu'nce  of  the  shock  is  small.  For  a  beam  with  irradiances  of  the  order 
of  a  few  kw/cm  the  vapor  Mach  number  is  subsonic  and  the  flow  field  beliaves 
as  a  quasi- steady,  subsonic  compressible  point  source.  The  pressure  at  the 
droplet  surface  is  related  to  the  vapor  Mach  number  iind  ambient  pressure  p^  by 

p  -  P  +  /2]  ,  (10) 

and  is  valid  for  M^<  1 . 0  and  R^^  CJ.  In  Eq.  (1  0),  Y  =  1  .  3  is  the  adi.ibatic  e xjionent 
of  water  vapor.  The  pressure  at  the  droplet  surface  relaxes  freim  the  initial  value 
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given  by  Eq.  (9)  to  a  value  slightly  less  than  one  atmosphere  in  a  few  acoustic 
times.  For  a  1  urn  droplet  this  occurs  in  times  less  than  0.  1  Us.  For  a  sub¬ 
sonic  point  source  the  velocity  decreases  and  the  pressure  increases  as  the 
distance  from  the  droplet  increases.  The  ambient  or  stagnation  pressure  p 

CD 

is  attained  when  the  radius  expands  to  twelve  times  the  initial  droplet  size. 

Once  ambient  pressure  is  reached  in  the  vapor  sphere  further  change  occurs 
by  diffusion  and  it  is  during  the  cooling  down  of  the  vapor  sphere  that  the  energy, 
the  sum  of  the  vapor  enthalpy  and  the  kinetic  energy  terms  in  Eq.  (3),  gets 
deposi-ed  into  the  atmosphere  as  heat.  The  droplet  temperature  and  mass  flux 
are  determined  as  a  function  of  Mach  number  from  Eq.  (1),  after  replacing  rn 
by  pv/RyTp  and  using  Eq.  (10)  for  the  pressure.  The  energy  balance  at  the  droplet 
surface,  Eq.  (3),  then  determines  the  absorbed  -rradiance  as  a  function  of  the 
mu-ss  flux. 

In  Fig.  3  we  present  the  change  in  droplet  temperature  as  a  function  of  the 

absorbed  irradiance  for  three  representative  droplets  and  two  values  of  the 

evaporation  coefficient.  The  droplets  were  initially  at  Z88*^K.  The  results  were 

generated  from  both  the  diffusion  and  the  conv'ective  analysis.  For  droplets 

heated  volumetrically  3Q^/4a  reduces  to  the  bulk  absorption  coi'fficient  Clp  .  At 

2 

3.  8  pm.'Tp  is  approximately  100/cm  and  the  irradiance  varies  from  10  w/cm 

b 

to  10  w/cm  in  Fig.  3.  At  10.6  pm,  Op  is  approximately  lOOC/cm  and  the  irradiance 
1  5  , 

varies  from  10  w/cm  to  10  w/cm.  The  solid  lines  arc  for  a  evaporation  coef¬ 
ficient  of  0.04  and  the  dashed  lines  for  an  evaporation  coefficient  of  1.0.  The 
10  n.icron  droplet,  w'ith  C  -  1.0,  is  the  only  case  that  exhibits  qua  si  -  cfiuilib  r  ium  • 
behaviour  over  a, wide  range  of  absorbed  i r radia nee s . 
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An  even  more  important  quantity  required  in  high-energy  laser  beam 
propagation  is  the  fraction  of  absorbed  energy  . 

f  =  [•O'r  +  m(Cp(Tp-TJ  +  v2/2)]/[kT^  +  m(L  +  Cp(Tp- T„)  +  v^/2)3  ,  (11) 

that  heats  the  atmosphere.  This  is  presented  in  Fig.  4  as  a  function  of  the 

absorbed  irradiance.  The  various  regimes  of  droplet  vaporization  are  clearly 

4  3 

illustrated  for  the  lO^m  droplet.  For  absorbed  irradiances  less  than  10  W /cm 
diffusive  transport  dominates.  In  this  regime  the  heat  conduction  term  KT^  is 
the  primary  contributor  to  f.  As  the  irradiance  is  increased  f  decreases,  at 
first,  because  of  the  exponential  decay  of  KT^  on  m.  The  relatively  flat  minimum 
represents  the  transition  between  diffusion  and  convection.  In  this  regime  the 

V 

vapor  enthalpy  term,  rn  Cp(Tp- T^^),  is  the  important  term  in  the  numerator  of 
Eq.  (11).  To  the  right  of  the  minimum  convection  dominates  and  f  increases 
because  the  stagnation  enthalpy,  the  sum  of  the  enthalpy  and  kinetic  energy 
terms  in  Eq.  (1  1),  is  increasing.  Droplets  other  than  lOpm  also  exhibit  the  same 
general  bchav'iour.  For  those  droplets  larger  than  10[am  the  miniinum  is  reached 
at  lower  absorbed  irradiances  and  for  those  droplets  smaller  than  10pm  the 
minimum  is  reached  at  higher  irradiances.  We  also  note  in  Fig.  4  the  difference 
in  results  for  the  two  evaporation  coefficients.  When  the  evaporation  coefficient 
is  1.0  at  least  50%  of  the  absorbed  energy  is  available  for  vaporizing  the  droplets 
When  tlie  evaporation  is  0.04  only  the  larger  droplets  arc  readily  vaporized. 

For  a  laser  propagating  in  the  atmosphere  the  appropriate  evaporation  coefficient 
is  0.04  and  tliese  results  indicate  that  water  droplets  have  longer  lifetimes 
than  prcvt'-  asly  expected. 


Mass  conservation  at  the  droplet  surfai^e. 


a  =  -(l  f(a.  I))QJ/4ppL,  (12) 

determines  the  change  in  radius.  Eq.  (12)  must  be  integrated  numerically 

because  of  the  complicated  dependence  of  f  on  o  and  I.  However  the  essential 

details  of  the  numerical  solutions  can  be  determined  from  the  behaviour  of 

droplets  volumetrically  heated  a  few  tens  of  degrees  above  the  ambient  tem- 

4  2 

perature.  For  an  evaporation  coefficient  of  0.04  and  3IQ^/4o<  10  w/cm  we  have 


1-f  =  1/(1.  56  +  2.  15/a{iJm)). 


(13) 


After  setting  Q^  =  4Gpa/3,  substituting  Eq.  (13)  into  Eq.(12)  and  integrating  we 
obtain  the  droplet  radius  in  microns  from 


1.  56  ln(ao/0)  +  2.  15(1/0  -  l/c^)  =  Iapt/3ppL  =  t.  (14) 

Results  for  initial  droplets  O^  =  0.  1,  1.0  and  10  microns  are  presented  in  Fig.  5 

as  a  function  of  the  normalized  time  T  =  Iapt/3PpL.  The  solid  lines  are  for  an 

4  3 

absor.jcd  irradiance  of  10  w/cm  and  the  dashed  lines  are  for  an  absorbed 
^  3 

irradiance  of  10  w/cm  .  The  actual  physical  time  for  the  curves  differ  by  a  * 

factor  of  a  hundred.  The  results  for  1  -  lO^w/cm^  were  determined  from 

Eq.{14)  and  numerical  integration  of  Eq.  (12)  was  used  to  determine  the  results 
6  3 

at  Icip  =  10  w/cm  .  The  only  significant  difference  in  the  results  for  the  two 
absorbed  irradiances  occurs  for  droplets  greater  than  Ipm.  Even  after  the  time 


12 


has  been  normalized  by  the  irradiance  a  10  )-im  droplet  initially  vaporizes  faster 
at  the  higher  iriadiances  because  more  of  the  energy  couples  into  vaporization. 
As  the  droplet  vaporizes  to  submicron  size  the  fraction  of  energy  coupled  into 
vaporization  becomes  equal  for  the  two  irradiances.  In  this  regime,  large  T  in 
in  Fig.  5,  all  droplets  are  approximately  equal  and  independent  of  the  initial  size 
and  absorbed  irradiance.  The  final  droplet  size  depends  only  on  the  normalized 
time  and  is  given  by  Eq.  (14)  . 
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IV.  PROPAGATION  ANALYSIS 


We  will  now  incorporate  the  single  droplet  results  into  an  analysis  for 
determining  the  irradiance  distribution  of  a  laser  beam  propagating  in  a 
hygroscopic  aerosol.  The  equation  describing  the  complex  ampHtutide  A  is 
the  paraxial  approximation  to  the  scalar  wave  equation 


2ikA^  +  V^A  +  2k^{TT^-l)pA  =  -ikf  )A  , 


(15) 


where  k  =  2tt/X  is  the  wave  number,  is  the  refractive  index  of  air,  is 

the  molecular  absorption  coefficient,  and  is  the  aerosol  extinction  coefficient. 

The  density,  p,  induced  by  molecular  and  aerosol  absorption  can  be  determined 

18 

from  a  simplified  version  of  the  hydrodynamics  of  heterogeneous  media  .  The 
variables  that  appear  in  the  continuity,  momentum  and  energy  equations  can 
be  regarded  as  averaged  over  a  gas  volume  containing  many  droplets.  Use 
of  the  average  values  requires  the  regions  of  immediate  influence  of  the  drop¬ 
lets  to  dissipate  rapidly  over  the  gas  volume.  The  vaporizing  droplets  act  as 
distributed  heat  and  mass  sources.  Multiplying  the  droplet  energy  balance, 

Eq.  (3)  by  4tto ^(0^)  and  integrating  over  Oq  gives  the  total  heat  added  to  the 
atmosphere.  Under  the  assumption  that  the  heating  occurs  at  constant  pressure 
the  steady  state  density  perturbation  induced  by  a  moving  aerosol  is 


Va,Pa 


V' t')dt'  ;u^  x-U(t-t') 


(16)  . 


where  Y  is  the  ideal  gas  constant  for  air,  p  is  the  ambient  pressure,  U  is 
00  oo 
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the  velocity  component  transverse  to  the  beam  and  is  the  heat  added 
to  the  atmosphere  by  the  aerosol.  The  total  extinction  coefficient  and  the 
heat  added  to  the  atmosphere  are  defined  by 

CD 

=  f  (17) 

o 

O) 

qg  =  y  na^(Oo.  I)f{o,I)Q^(a(o^.I))n(o^)dOQ 
o 

where  n(  a^)  is  the  initial  droplet  distribution  function.  The  Mie  absorption 

19 

coefficient  can  be  calculated  using  the  following  analytical  approximation 

Qa  =  l-(l-(l  +  2apa)exp(-2Q^rr))/2(apCT)^;  Op  =  2kri.  (19) 

However  the  extinction  coefficient  should  be  calculated  using  the  exact  theory 

®  2 

Qg  =  X/ (4n+2)Re(a|^  +  bj,)/ag  ;  =  Zna/X  (20) 

®  0 

Evaluation  of  the  irradiance  dependence  of  and  requires  solution  of 
E  q.  (12)  for  those  droplets  in  the  size  range  specified  by  n(a  ). 

This  completes  the  formulation  for  aerosol  thermal  blooming.  The  case 
of  non- vaporizing  particles  such  as  carbon  particles  is  obtained  from  the 
above  analysis  by  setting  f(o,  I)=l  and  o=  in  Eqs.  (17)  and  (18).  The  results 
for  the  extinction  coefficient  and  the  heat  addition  term  q  are  sensitive  to 
the  laser  wavelength  and  size  distribution  of  the  droph-ts. 
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Two  dropU't  distributicMi  functions  that  art-  much  usi-d  fur  modiMlinp 


aerosols  are  the  modified  gamma  distribution  of  Deirmendjian  and  the 

1  b 

power  law  distribution  of  Junge 


a  b  -  a 

n(a^)  =  N(Oo/a^)  exp(-b(OQ/o^)  )  ;  n(ao)  = 


(21) 


The  parameters,  in  the  modified  gamma  distribution  for  various  types  of 
aerosols,  are 


haze 

a=  2 

b  =  4 

c  =  .  5 

0  =  .  07  (-im 

m 

Oj^  =  .  07pm 

rain 

a=  2 

b  =  4 

c=  .  5 

clouds 

II 

00 

b  =  6 

c  =  1 

0  =4.0  pm 

m 


The  parameters  in  the  power  law  distribution  for  various  types  of  aerosols 
are 

haze  3<a<4  .01  <0^'^  10 

fogs-smokes  a=2  ,l<O|^<20 

dust  clouds  a=2  .4<Oq<’4 


To  complete  the  analysis  we  need  only  specify  the  real  and  imaginary 
parts  of  the  index  of  refraction  and  typical  mass  extinction  coefficients  for 
battlefield  aerosols.  These  aerosols  include  smokes  that  are  hygroscopic, 
dust  clouds,  carbon  smokes  such  as  diesel  exhaust  and  fog  oil.  At  10.  6  pm 
these  properties  are  summarized  as  follows: 

Red  Phosphorus 
fog  Oil 
Carbon  Smokes 
Dust  Clouds 


nr 

ni 

a^[m 

1.8 

0.  389 

0.  244 

0.  218 

1. 48 

0.  163 

0.  168 

0.  167 

3.  3 

1. 70 

0.  80 

>.  9ap 

0.  20 

>. 
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For  non- vaporizing  aerosols  these  results  can  be  incorporated  into  the 
Army  propagation  codes  by  adding  q^,  with  f=  1  and  a  =  cIq  ,  to  the  thermal 
blooming  contribution  and  adding  with  a=  Oq,  to  the  exponential  term 
(i.  e.  exp(- (a^  +  a^)z)  ).  For  vaporizing  aerosols  this  is  still  correct  except 
now  we  must  compute  f  and  a  with  an  aerosol  subroutine  that  is  called  for 
each  value  of  x,  y,  and  time  t.  This  subroutine  then  computes  0.^  and 
using  the  local  values  of  the  irradiance  I(x,  y,  z),  0(0^,  t),  and  f(r,  I),  In 
order  that  this  maybe  accomplished  the  change  in  droplet  size  for  each 
droplet  in  the  droplet  distribution  function  must  be  stored  for  each  x  and  y. 

The  Mie  absorption  and  extinction  using  the  new,  vaporized  ,  droplet  sizes 
is  then  used  to  compute  and  qg.  This  value  is  then  passed  to  the  main 
propagation  program. 

For  pulsed  lasers  this  indicates  that  any  relevant  propagation  calculation 
must  be  four  dimensional.  As  an  example  of  the  magnitude  of  the  calculation 
we  must  compute  and  store  the  change  in  droplet  size  for  about  200  droplets 
at  each  value  of  x,  y,  and  t  ,  determine  the  200  Mie  extinction  coefficients 
for  these  vaporized  droplets,  and  do  this  for  each  transverse  mesh  point 
in  the  propagation  calculation.  For  a  transverse  mesh  point  grid  64x64  and  20 
z-steps  about  10  million  Mie  extinction  calculations  are  required  in  the  typical 
propagation  calculation.  Since  Mie  scattering  calculations  are  well  known  to  be 
time  consuming  development  of  efficient  algoritms  are  mandatory.  The  only 
other  approach  is  to  develop  rational  (analytical  expressions)  approximations 
to  the  Mie  coefficients.  Presently,  accurate  expressions  over  the  entire  droplet 
size  range  are  not  known.  Until  this  is  accomplised  full  numerical  treatment 
of  nonlinear  aerosol  thermal  blooming  is  not  possible. 
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Mv=  v/(>'vRvTp)'/2 

Fig.l  Macii  number  ciejjendcnce  of  the  ratio  of  the* 
saturated  vapor  pressure  to  tiic  actual  vapor  pressure 
Amount  of  superheat  is  limited  by  the  tensile  strength 
of  tlie  droplet.  Also  shown  are  the  asymptotes  where 
collisions  of  the  vapor  molecules  with  the  droplet  are 
small. 


ions  for  the  diffusion 
Is.  Int('rfacial  discijn- 
order  of  a  mean  free 


3IQa/4c7  (W/cm3) 

Fig.  3  Dependence  of  the  peak  droplet  temperature 
on  the  absorbed  irradiance.  For  droplets  that 
absorb  volumetrically  the  abscissa  becomes  I&_. 


31Qa/4cr  (W/cm^) 


Fi^.  4  Fraction  of  the  absorbed  energy  that  is 
deposited  into  the  atmosphere  as  heat.  The  re¬ 
mainder  of  the  absorbed  energy  vaporizes  the 
droplet. 


\ 
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DROPLET 


I'iu.  5  Droplet  lifetimes 
i  r  r  fl  ia  n  c  e  .  Also  sliowr 
life  limes  unrler  tlie  a  dir 


